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Abstract 

Recently, gravitational microlensing has been investigated in the framework of the weak field 
limit of fourth order gravity theory. However, solar system data (i.e. planetary periods and light 
bending) can be used to put strong constraints on the parameters of this class of gravity theories. 
We find that these parameters must be very close to those corresponding to the Newtonian limit 
of the theory. 
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I. INTRODUCTION 



Since many years different alternative approaches to gravity have been proposed in the 



literature such as MOND 



theories 



1 fl 



sj3,y,3, 



2j, scalar-tensor 



3| , conformal 4] , Yukawa-like corrected gravity 



and so on (see papers [8| for reviews). Very recently, it has been proposed 
111 ], in the framework of higher order theories of gravity - also referred to as f(R) 



theories - a modification of the gravity action with the form 



A 



d 4 x^{f(R)+£ rt 



(1) 



where f(R) is a generic function of the Ricci scalar curvature and C m is the standard matter 
Lagrangian. For example, if f(R) = R + 2 A the theory coincides with General Relativity 



(GR) with the A term. In particular, Capozziello et al. 1CJ, 111] considered power law 



function f(R) theories of the form f(R) = f()R n - As a result, in the weak field limit 13], 



the gravitational potential is found to be lOj, llll ]: 



$0) 



Gm 
"27 



where 



P 



12n 2 -In - I- V36n 4 + 12n 3 - 83n 2 + 50n + 1 



(2) 



(3) 



Qn 2 - An + 2 

The dependence of the (3 parameter on the n power is shown in Fig. [TJ Of course, for 
n — > oo it follows /3 — *■ 1, while for n — 1 the parameter (3 reduces to zero and the Newtonian 
gravitational field is recovered. On the other hand, while (3 is a universal parameter, r c in 
principle is an arbitrary parameter, depending on the considered system and its typical scale. 
Consider for example the Sun as the source of the gravitational field and the Earth as the 
test particle. Since Earth velocity is ~ 30 km s _1 , it has been found that the parameter 
r c varies in the range ~ 1 — 10 4 AU. Once r c and (3 has been fixed, Capozziello et al. [h]] 
used them to study deviations from the standard Paczynski light curve f or g ravitational 
micro lensing [3] and claimed that the implied deviation can be measured 15|. It is clear 



that for gravitational microlensing one could detect observational differences between GR 
and an alternative theory (the fourth order gravity in particular), so that one should have 
different potentials at the scale Re (the Einstein radius) of the gravitational microlensing 
. For the Galactic microlensing case Re is about 1 AU. This is a reason why the authors 
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FIG. 1: The parameter (5 as a function of n for fourth order gravity. 



|10l ] have selected r c at a level of astronomical units to obtain observable signatures for non- 
vanishing (3. The aim of the present paper is to show that solar system data (light bending 
and planetary periods) put extremely strong constraints on both r c and f3 parameters making 
this alternative theory of gravity not so attracting. 

II. SOLAR SYSTEM CONSTRAINTS 
A. Light Bending Constraints 

As stated above, we now discuss some observational consequences of the fourth order 
gravity, depending on the choice of the parameters (3 and r c . 

A constraint on the proposed theory can be derived by considering the light bending effect 
in the Sun limb. It is well-known that in the parameterized post-Newtonian formalism the 
bending angle through which a electromagnetic light ray from a distant source is deflected 
by a body with mass m is jsl, Q 



where b is the impact parameter, cb is the solar elongation angle (between the Sun and the 




(4) 



3 



source as viewed from Earth) and 7 is the post-Newtonian parameter. For GR, 7 = 1 and 
for light rays at Sun's limb, #gr = 1.75". Recently, Shapiro et al. measured the bending 
angles for distant compact sources and concluded that light bending angles follow GR with 
a very high precision (7 = 0.9998 ± 0.0004). In other words, this means that the deflection 
of the light path is well described by the GR theory. In particular, as radio observations 
of distant sources have shown [18], the observed and expected bending angles are related 
by 9 bs = (1.001 ± 0.001)6'gr. In the framework of the fourth order gravity theory, the 
deflection angle of light rays at Sun's limb depends on both the parameters f3 and r c . We 
explore this dependence in Fig. [2, by requiring that the expected value for the bending 
angle is, at least, within 2a (grey region) or within 5a (light grey region) the observed 
value. Inspecting the same figure, it is clear that only /5-values nearby zero (corresponding 
to a completely arbitrariness of r c ) are consistent with the observed deflection angles. We 
therefore emphasize that /3-values considered in [k]] and 
are ruled out by light deflection data. 



ll| (i.e. = 0.25,0.43,0.58,0.75) 



B. Planetary Constraints 



A stronger constraint on the fourth order gravity theory can be obtained from the motion 
of the solar system planets. Let us consider as a toy model a planet moving on circular orbit 
(of radius r) around the Sun. From Eq. (TSJ), the planet acceleration a = — <9<3>(r)/<9r is given 
by 

Accordingly, the planetary circular velocity v can be evaluated and, in turn, the orbital 
period P is given by 

a ai -1/2 



P = P K V2 



where Pk = [47r 2 r 3 / (Gm)} 1 / 2 is the usual Keplerian period. In order to compare the orbital 
period predicted by the fourth order theory with the Solar System observations, let us define 
the quantity 

AP I p — p„\ 

l/(Ar c )-l| (7) 
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FIG. 2: Constraints on the fourth order theory parameters {(5 and r c ) arising from the deflection 
angle of light rays close to the solar limb. The grey and light grey regions embed the part of the 
parameter space allowed by solar system observations at the 2a and 5a confidence level, respec- 
tively. It is noticing that for scale reason we have not plotted values of r c up to 10 4 AU. For these 
cases, the observations can be restored only for (3 — * 0. 

being r c ) the factor appearing on the right hand side of equation (jSJ) and multiplying 
the usual Keplerian period. 

There is a question about a possibility to satisfy the planetary period condition - van- 
ishing the Eq. (JTj) - with f3 parameter which is significantly different from zero. Vanishing 
the right hand side of Eq. (J7]) we obtain the relation 



so that Eq. (jSJ) should be satisfied for all the planetary radii. This is obviously impossible 
since the fourth order theory defines (3 as a parameter, while the specific system under 
consideration (the Solar system in our case) allows us to specify the r c parameter. Hence 
the right hand side of Eq. (jSJ) is fixed for the Solar system, implying that it is impossible to 
satisfy Eq. (jSJ) even with two (or more) different planetary radii. 

Just for illustration we present the function f(/3, dependence on (3 parameter for 

fixed r c and planetary radii r (see Fig. [3]). As one can see from Eq. (jSJ) (and Fig. [3] as 



In r = In r r — 




(8) 
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FIG. 3: The orbital period in units of the Keplerian one (the function f(f3,r c )) is given, as a 
function of ft, in the case of Mercury (dashed line), Venus (dotted line), Earth (solid line) for 
different r c (it is clear that f(ft,r c ) — * V% for ft — ► 1). If r > r c there is ft G (0,1) satisfying 
Eq. (jSJ , but the ft values depend on fixed r c and r and they are different for a fixed r c and different 
r, so that it is impossible to satisfy Eq. (jHJ) if the number of planets is more than one. Moreover, 
if we have at least one radius r < r c , there is no solution of Eq. (jHJ). 

well) for each planetary radius r > r c there is ft G (0, 1) satisfying Eq. (jHJ), but the ft value 
depends on fixed r c and r, so that they are different for a fixed r c and different r. Moreover, 
if we have at least one radius r < r c , there is no solution of Eq. (jHJ). Both cases imply that 
the ft parameter should be around zero. 

In Fig. HJ (left panel), the factor f(ft,r c ) is given as a function of ft for the two limiting 
values of r c , 1 AU (dashed line) and ~ 10 4 AU (solid line), considered in [loj]. As one can 
note, only for ft approaching zero it is expected to recover the value of the Keplerian period. 
In the above mentioned figure, the calculation has been performed for the Earth orbit (i.e. 
r = 1 AU). 

Current observations allow also to evaluate the distances between the Sun and the planets 
of the Solar System with a great accuracy. In particular, differences in the heliocentric 
distances do not exceed 10 km for Jupiter and amount to 180, 410, 1200 and 14000 km 



for Saturn, Uranus, Neptune and Pluto, respectively |19|]. Errors in the semi- major axes 
of the inner planets are even smaller (see e.g. Table 2 in 20]) so that the relative error 
in the orbital period determination is extremely low. As an example, the orbital period of 
Earth is T = 365.256363051 days with an error of AT = 5.0 x 10 -10 days, corresponding 
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FIG. 4: The factor /(/?, r c ) is given as a function of (3 (left panel) for the two limiting values of r c , 1 
AU (dashed line) and ~ 10 4 AU (solid line), respectively. As one can note, only for j3 approaching 
it is expected to recover the value of the Keplerian period. Here, the calculation have been 
performed at Earth (i.e. r = 1 AU). In the middle and right panel, the quantity /(/3, r c ) — 1 is 
given as a function of (3 for r c = 10 AU and r c = 1 AU. Note that only for values of (3 close to 
the Solar System observation can be restored (see text for more details). 

to a relative error of AT/T less than 10~ 12 . These values can be used in order to constrain 
the possible values of both the parameters (3 and r c introduced by the fourth order gravity 
theory. This can be done by requiring that AP/Px < AT/T so that, in the case of Earth, 
\f{/3, r c ) — 1| < 10~ 12 which can be solved with respect to (3 once the r c parameter has been 
fixed to some value. For r c = 1 AU and r c = 10 4 AU (i.e. the two limiting cases considered by 
Capozziello et al. [loj]) we find the allowed upper limits on the f3 parameter to be 4.0 x 10~ 12 
and 3.9 x 10~ 13 , respectively (since AP/Pk = Af3[— 1 + In (r/r c )]/4). These results can also 
be inferred from the middle and right panels of Fig. |H 

A more precise analysis which takes into account the planetary semi-major axes and 
eccentricities leads to variations of at most a few percent with respect to the results in 
Fig. HI since the planetary orbits are nearly circular. Therefore, in spite of the fact that 
orbital periods of planets are not generally used to test alternative theories of gravity (since 
it is taken for granted that the weak field approximation of these theories gives the Newtonian 
limit), we found that these data are important to constrain parameters of the fourth order 
gravity theory. 
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III. DISCUSSION 



GR and Newtonian theory (as its weak field limit) were verified by a very precise way 
at different scales. There are observational data which constrain parameters of alternative 
theories as well. As a result, the parameter (3 of fourth order gravity should be very close to 
zero (it means that the gravitational theory should be very close to GR). In particular, the (5 



parameter values considered for microlensing [10], for rotation curves 



ll| and cosmological 



SN type la 21[ are ruled out by solar system data. 



No doubt that one could also derive further constraints on the fourth order gravity theory 
by analy zing other physical phenomena such as Shapiro time delay, frequency shift of radio 



photons [22j, laser ranging for distant objects in the solar system, deviations of trajectories 
of celestial bodies from ellipses, parabolas and hyperbolas and so on. But our aim was only 
to show that only f3 ~ values are not in contradiction with solar system data in spite of the 
fact that there are a lot of speculations to fit observational data with /3 values significantly 
different from zero. 
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